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, What is claimed is: 

\i/ \ j Jp^ 1. An apparatus for computing multiple integral of a multidimensional integrand 
function A tcr^e integrated with using a vector map f with unbounded support which 
converts an m (r^>l) -dimensional vector having real number components into an m- 
5 dimensional vectonhaving real number components, by which a multidimensional density 
function p for th^Iimk^ distribution resulting from repeatedly applying the map f to an 
m-dimensional vector ims analytically solvable, said apparatus comprising: 

a first storage unit which stores an m-dimensional vector x = (Xi, . . . , x^J; 
a second storage unit which stores a scalar value w; 
10 a first computing unit wnich computes a vector x' = f(x) = (x\, x'25 . . . r^'m) resulting 
from applying said vector map Ato said vector x being stored in said first storage unit; 

a second computing unit which computes a scalar value w' = A(x)/p(x) based on 
said vector x being stored in said fir^t storage unit and said scalar value w being stored in 
said second storage unit; 
15 an update unit which updates the ^ue stored in said first storage unit by storing 
said vector x' computed by said first commiting unit on said first storage unit, and updates 
the value in said second storage unit by adding said scalar value w' computed by said 
second computing unit to a value to be storedVn said second storage unit; and 

an output unit which computes a scalar vame s = w/(c+l) based on said scalar value 
20 w being stored in said second storage unit when tke number of update times by said 
update unit becomes c (c>l), and outputs said scalsu^ value s as a result of the multiple 
integral. 

2. The apparatus according to claim 1, wherelh said scalar value stored in said 
second storage unit first is a result from dividing a value Vesulting from applying said 
25 function A to said m-dimensional vector stored in said firsKstorage unit first by a value 
resulting from applying said density function p to said m-dii\ensional vector stored in 
said first storage unit first. 
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3. nrhe apparatus according to claim 1, wherein said scalar value stored in said 
second storage unit first is 0, and 

said output unit computes a scalar value s' = w/c, and outputs said scalar value s' as 
the result of the multiple integral instead of said scalar value s. 
5 4. The apparatus according to claim 1 further comprising: 

a convergence \ate obtainer which obtains convergence rate of scalar values 
sequentially output by^aid output unit while varying said number of update times c for 
each of plural vector m^s gi, g25 . . . , g^ (k>2) which are prepared as said vector map f; 
a vector map select^- which refers to the convergence rates obtained by said 
1 0 convergence rate obtainer, £^d selects a vector map g^ (l<h<k) which shows fastest 
convergence rate; and 

an output controller whicflk controls said output unit to output said scalar value with 
using said vector map gi, as said vector map f and the number of update times c' (c'>c) 
instead of said number of update times c. 
15 5, The apparatus according do claim 1, wherein a multidimensional density 
function p representing the limiting distmbution of a vector sequence 

u, f(u), f(f(^), f(f(f(u)))), . . . 
resulting from applying said vector map f tt\a predetermined m-dimensional vector u = 
(u,, Uj, . . . , lO for equal to or greater than 0 nines, satisfies the following property of: 

20 p(u) - ni=oV(Ui); 

PiCuO-cJur'""^ toru.--a): 
p^CuO-cJur'^") forV^+co: 
(a>0, l<i<m, c_>0, c>0) 
6. The apparatus according to claim 5, whereii^aid vector map f is defined as 

25 f(u) = (f,(uO, f2(uj, . . . , 4(1 

by a function fi(t) = gi(di tyd^ (di>0) which is defined \n l<i<m, and said map gi is 
defined by any one of the following maps cpj (l<j<8) and a natural number n^ (n>2), as 
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follows: 

gi((Pj(e)) = cpj(nie); 
(PiV) = -sgn(tan0) / |tan0|^'^ 
<p2(6)\-sgn(tan0) x ItanGp^^; 
93(9) = -\gn(cose) / Itaner^ 
(P4(0) - -s^/cosB) X ItanGp'^ 
(P5(0) = sgn(cds0)/|tan0r^ 
(PeCG) = sgn(cos0^ X itanGr^ 
(P7(0) = sgn(sin0) /^an0|^^"; 
(PsCG) = sgn(sin0) X ithnGp^- 
sgn(t) =1 for t > 0; 
sgn(t) = 0 for t = 0; 
sgn(t) = -l fort<0 

7. The apparatus accordin^o claim 5 further comprising: 
a convergence rate obtainer whicla defines said map f for each of plural positive 

numbers (h, . , . , (k^2) prepared a^an invariable a, and obtains convergence rates of 
the scalar values sequentially output by sai\j output unit while varying said number of 
update times c; 

a positive number selector which refers tlk convergence rates obtained by said 
convergence rate obtainer, and selects an integer \jh (l<h<k) which shows the fastest 
convergence rate; and 

an output controller which defines said map f v)ith using said positive number as 
said invariable a, and controls said output unit to output said scalar values with using the 
number of update times c' (c*>c) instead of said number ©f update times c. 

8. The apparatus according to claim 6 further comprising: 
a convergence rate obtainer which defines said map g^ with using plural ones of said 

maps (pj, and obtains convergence rates of the scalar values sequentially output by said 
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output unit while varying said number of update times c; 

\ a map selector which refers to the convergence rates obtained by said convergence 
rate ootainer, and selects one of said maps (pj which shows the fastest convergence rate; 
and \ 

5 an outpqt controller which defines said map gi with using said map (pj selected by 
said map selectoi^ and controls said output unit to output the scalar values with using the 
number of update iknes c' (c'>c) instead of said number of update times c. 

9, • The apparatus according to claim 6 further comprising: 

a convergence ratei obtainer which defines said map g^ relating to each of plural 
1 0 natural numbers Pi, P2, . . V Pk (k^2) as said natural numbers n^, and obtains convergence 
rates of the scalar values secmentially output by said output unit while varying said 
number of update times c; \ 

a natural number selector wnich refers to the convergence rates obtained by said 
convergence rate obtainer, and selects a natural number p^ (l<h<k) which shows the 
1 5 fastest convergence rate; and \ 

an output controller which defines said natural number map g^ with using said 

natural number p^ as said natural number n\ and controls said output unit to output the 

scalar values with using the number of updatestimes c' (c'>c) instead of said number of 

update times c. V 

20 10. The apparatus according to claim 1, v^erein said output unit computes said 

scalar value s each time said update unit carries out update, compares said latest scalar 

value s with the former scalar value which is computed at former update, and outputs said 

latest scalar value s if a result of the comparison satisfies^a predetermined condition for 

terminating the computation. 

25 11. A method for computing multiple integral of a multidimensional integrand 

\ 

function A to be integrated with using a vector map f with unbounded support which 
converts an m (m>l) -dimensional vector having real number components into an m- 



31 



dimensional vector having real number components, by w^hich a multidimensional density 
function p for the limiting distribution resulting from repeatedly applying the map f to an 
m-dirrwisional vector u is analytically solvable, 

a first storage unit w^hich stores an m-dimensional vector x = (x^, X2, . . . , xj, and 
5 a second storage unit which stores a scalar value w, 
said method comprising the steps of: 

computing aWector x' = f(x) = (x\, x'2, . . . , x'^J resulting from applying said vector 
map f to said vectorV being stored in said first storage unit; 

computing a scalar value w' = A(x)/p(x) based on said vector x being stored in said 
1 0 first storage unit and said scalar value w being stored in said second storage unit; 

updating the value stored in said first storage unit by storing said vector x* computed 
by^aid first computing unit on said first storage unit, and updating the value in said 
second storage unit by adding s^d scalar value w' computed b>\^id^ second computing 
unit to a value to be stored on saio^second storage unit; and 
1 5 computing a scalar value s = v^(c-M) based on said scalar value w being stored in 
said second storage unit when the number of update times by said update unit becomes c 
(c>l), and outputting said scalar value sVs a result of the multiple integral. 

12. The method according to claim 1 1, wherein said scalar value stored in said 
second storage unit first is a result from divimng a value resulting from applying said 

20 function A to said m-dimensional vector storea in said first storage unit first by a value 
resulting from applying said density function p to^said m-dimensional vector stored in 
said first storage unit first. 

13. The method according to claim 11, wherd^n said scalar value stored in said 
second storage unit first is 0, and 

25 said outputting step computes a scalar value s' = w/c\ and outputs said scalar value s' 
as the result of the multiple integral instead of said scalar value s. 

14. The method according to claim 1 1 further comprising the steps of: 




32 

olkaining convergence rate of scalar values sequentially output by said outputting 
step whileWarying said number of update times c for each of plural vector maps g^, 
§2. • ' . gk which are prepared as said vector map f; and 

referring\p the convergence rates obtained by said convergence rate obtainer, and 
5 selecting a vecto\map (l<h<k) which shows fastest convergence rate, and 

said outputting step outputs said scalar value with using said vector map g^ as said 
vector map f and the number of update times c' (c'>c) instead of said number of update 
times c. 

15. The method According to claim 1 1, wherein a multidimensional density 
1 0 function p representing the Ijmiting distribution of a vector sequence 

.u,f(uXf(f(u)Xf(f(f(u))X... 
resulting from applying said vector map f to a predetermined m-dimensional vector u = 
(ui, U2, . . . , Mra) for equal to or greater than 0 times, satisfies the following property of: 
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P.(^)-c_>lur-> for^---a.; 
Pi(Ui) - |u^^'^^> for Ui ^ +0); 
(a>0, l<i<V c_i>0, c,i>0) 
16. The method according to claim 13,^ wherein said vector map f is defined as 

f(u) = (f,(uO,f,(ujV..,4,(uJ) 

by a function fi(t) = gi(di t)/di (di>0) which is defined in l<i<m, and said map g; is 
defined by any one of the following maps (pj (l<j<8\and a natural number n^ (n>2), as 
follows; 

gi((Pj(e)) = (Pj(nie); 

9,(9) = -sgn(tane) / Itanei""; 

92(8) = -sgn(tane) x Itanei""; 

93(6) = -sgn(cose) / Itanei""; 

94(e) = -sgn(cose) X Itanei""; 
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cpsXe) = sgn(cos0) / Itaner^; 
(P6(e\ = sgn(cose) X ItanGI^'^ 
(pXe) V sgn(sine) / Itaner^; 
(9) =Vgn(sine) x |tan0r^ 
5 sgn(t) = iv fort>0; 
sgn(t) = 0 \ for t = 0; 
sgn(t)=-l fort<0 

17. The memod according to claim 15 further comprising the step of: 

defining said map f for each of plural positive numbers Qi, • • • , qk (k>2) prepared 
1 0 as an invariable a, andpbtaining convergence rates of the scalar values sequentially 
output by said output unVt while varying said number of update times c; 

referring to the obtamed convergence rates, and selecting a positive number q^, 
(l<h<k) which shows the fastest convergence rate; and 

defining said map f witA using said positive number q^, as said invariable a, and 
1 5 controlling output of said scalaK values with using the number of update times c' (c'>c) 
instead of said number of update \imes c. 

18. The method according \o claim 16 further comprising the steps of: 
defining said map g^ with usingplural ones of said maps (pj, and obtaining 

convergence rates of the scalar values sequentially output by said output step while 
20 varying said number of update times c; 

referring to the obtained convergenceVates, and selecting one of said maps cpj which 
shows the fastest convergence rate; and 

defining said map gi with using said seleoted map cpj selected, and controlling output 
of the scalar values with using the number of update times c' (c'>c) instead of said number 
25 of update times c. 

19. The method according to claim 16 further comprising the steps of: 
defining said map gj relating to each of plural natural numbers Pi, P2, • • , Pk (k>2) 



/ as said natural numbers n^, and obtaining convergence rates of the scalar values 

sequentially output by said outputting step while varying said number of update times c; 

\eferring to the obtained convergence rates, and selecting a natural number p^^ 
(l<h<k)Yhich shows the fastest convergence rate; and 
5 definmg said natural number map g^ with using said natural number p^, as said 
natural numbW n^, and controlling output of the scalar values with using the number of 
update times c'\c'>c) instead of said number of update times c, 

20. The method according to claim 11, wherein said outputting step computes 
said scalar value s eVch time said updating step carries out update, compares said latest 

1 0 scalar value s with theVormer scalar value which is computed at former update, and 
outputs said latest scalarValue s if a result of the comparison satisfies a predetermined 
condition for terminating tne computation. 

21. A computer readable recording medium storing a program for computing 
multiple integral of a multidimensional integrand function A to be integrated with using a 

15 vector map f with unbounded support which converts an m (m>l) -dimensional vector 
having real number components intoi an m-dimensional vector having real number 
components, by which a multidimensional density function p for the limiting distribution 
resulting from repeatedly applying the map f to an m-dimensional vector u is analytically 
solvable, said program causes a computer tsd function as: 

20 a first storage unit which stores an m-dVnensional vector x = (Xi, X2, . . . , x^; 
a second storage unit which stores a scalar value w; 

a first computing unit which computes a veVtor x' = f(x) = (x'l, x'2, ^ . . , x'^J resulting 
from applying said vector map f to said vector x beVg stored in said first storage unit; 

a second computing unit which computes a scalV value w' = A(x)/p(x) based on 
25 said vector x being stored in said first storage unit and sVd scalar value w being stored in 
said second storage unit; \ 

an update unit which updates the value stored in said fiVst storage unit by storing 
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Said vector x' computed by said first computing unit on said first storage unit, and updates 
the\vaiue in said second storage unit by adding said scalar value W computed by said 
second computing unit to a value to be stored on said second storage unit; and 

an\)utput unit which computes a scalar value s = w/(c+l) based on said scalar value 
w being stored in said second storage unit when the number of update times by said 
update unit becomes c (c>l), and outputs said scalar value s as a result of the multiple 
integral. 

22. The recording medium according to claim 21, wherein said scalar value 
stored in said seconovstorage unit first is a result from dividing a value resulting from 
applying said fiinctionV to said m-dimensional vector stored in said first storage unit first 
by a value resulting fromV)pIying said density function p to said m-dimensional vector 
stored in said first storage unit first. 

23. The recording medWm according to claim 21, wherein said scalar value 
stored in said second storage unitVrst is 0, and 

said output unit computes a scalar value s' = w/c, and outputs said scalar value s' as 
the result of the multiple integral insteVi of said scalar value s. 

24. The recording medium according to claim 21, wherein said program further 
causes said computer to function as: ' 

a convergence rate obtainer which obtaihs convergence rate of scalar values 
sequentially output by said output unit while varVing said number of update times c for 
each of plural vector maps gt, g2, . . . , g^ (k>2) whffeh are prepared as said vector map f; 

a vector map selector which refers to the convergence rates obtained by said 
convergence rate obtainer, and selects a vector map g^ (l^h<k:) which shows fastest 
convergence rate; and 

an output controller which controls said output unit to oii^ut said scalar value with 
using said vector map g^^ as said vector map f and the number of\pdate times c' (c'>c) 
instead of said number of update times c. 
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' ^y^/^ TT^e recording medium according to claim 21, wherein a multidimensional 

density function p representing the limiting distribution of a vector sequence 

u, f(u), f(f(u)), f(f(f(u))), ... 
resulting from applying said vector map f to a predetermined m-dimensional vector u = 
5 (u,, Uj, . . . , u^ for equal to or greater than 0 times, satisfies the following property of: 

P(u) = n^o" P,(Ui); 
p,(uD~c_Jur"^ for ^ ->-«); 
Pi(Ui) ~ c_i luj-^'^") for Ui ^ +co; 
(a>0, l<i<m, c_i>0, c^,>0) 
1 0 26. The recorojng medium according to claim 25, wherein said vector map f is 
defined as 

f(u) = (f,(u,), f,(uj, . . . , f;„(uj) 

by a function fi(t) = gM^ t)/di (di>0) which is defined in l<i<m, and said map g; is 
defined by any one of the following maps (pj (l<j<8) and a natural number (n>2), as 
15 follows: 

&((Pj(0)) = (p/niG); 
(p,(e) = -sgn(tane) / ItanOr^ 
(P2(e) = -sgn(tane) x |tan0|"^ 
93(9) = -sgn(cos9) / Itanei""; 
20 (p4(e) = -sgn(cose) X Itanei""; 
(P3(e) = sgn(cose)/|tane|"-; 
(PsCe) = sgn(cose) X Itanei""; 
(P7(e) = sgn(sin0) / ItanGI""; 
(P8(e) = sgn(sine) xltanGI'""; 
25 sgn(t) = 1 for t > 0; 
sgn(t) = 0 for t = 0; 
sgn(t)=-l fort<0 
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7\ The recording medium according to claim 25, wherein said program further 

(ry \ 

causes said computer to function as: 

a coiWergence rate obtainer which defines said map f for each of plural positive 
numbers qi,>g2, . - - , % (k>2) prepared as an invariable a, and obtains convergence rates of 
5 the scalar valiJes sequentially output by said output unit while varying said number of 
update times c; 

an integer selector which refers the convergence rates obtained by said convergence 
rate obtainer, and se^cts a positive number q^, (l<h<k) which shows the fastest 
convergence rate; and 

10 an output controller which defines said map f with using said positive number as 

\ 

said invariable a, and controls said output unit to output said scalar values with using the 



number of update times c' (^5>c) instead of said number of update times c. 

28. The recording medium according to claim 26, wherein said program further 
causes said computer to function bs: 

15 a convergence rate obtainer which defines said map gj with using plural ones of said 
maps cpj, and obtains convergence ratesspf the scalar values sequentially output by said 
output unit while varying said number oAipdate times c; 

a map selector which refers to the con\^rgence rates obtained by said convergence 
rate obtainer, and selects one of said maps (pj wiich shows the fastest convergence rate; 

20 and 

an output controller which defines said map A with using said map cpj selected by 
said map selector, and controls said output unit to ourput the scalar values with using the 
number of update times c' (c'>c) instead of said numbenof update times c. 

29. The recording medium according to claim 26,\vherein said program further 
25 causes said computer to function as: 

a convergence rate obtainer which defines said map g^ refating to each of plural 
natural numbers Pt, P25 - • * > Pk (k>2) as said natural numbers ni,\^id^^ convergence 
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rates of the scalar values sequentially output by said output unit while varying said 
nuhiber of update times c; 

aNimtural number selector which refers to the convergence rates obtained by said 
convergencb^rate obtainer, and selects a natural number p^ (l<h<k) which shows the 
fastest convergence rate; and 

an output controller which defines said natural number map g^ with using said 
natural number p^, as said natural number n^, and controls said output unit to output the 
scalar values with using me number of update times c' (c'>c) instead of said number of 
update times c. \ 

30. The recording medium according to claim 21, wherein said output unit 
computes said scalar value s each time said update unit carries out update, compares said 
latest scalar value s with the former scalW value which is computed at former update, and 
outputs said latest scalar value s if a result of the comparison satisfies a predetermined 
condition for terminating the computation. \ 

31. The recording medium according to cTaim 21, wherein said recording medium 
is a compact disc, a floppy disk, a hard disk, a magneto-optical disk, a digital versatile 
disc, a magnetic tape, or a semiconductor memory, \ 



